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Abstract. In the present paper, we estimate the convergence rate in the 
Janashia-Lagvilava spectral factorization algorithm (see Studia Mathematica, 
137, 1999, 93-100) under the restriction on a spectral density matrix that its 
inverse is integrable. 
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1. Introduction. Wiener's spectral factorization theorem [7], [8] (see 
also [2], [1]) asserts that a positive definite matrix- valued function 

(1) S(z) = (Mz))l J=1 , 

fij G Li(T), on the unit circle of the complex plane with integrable 
logarithm of the determinant 

(2) log det (S(z)) e L X (T), 
admits the representation 

(3) S(z)= X + (z)(x + (z)T, 
where 

oo 

(4) X + (*) = E^> 

k=0 

\z\ < 1, 7^ are rxr matrix coefficients, is an analytic matrix-function 
with entries from the Hardy space if 2, X + { z ) £ H2, and the determinant 
of which is an outer function. 

The relation (3) is assumed to hold a.e. on the unit circle T and 
— (x + ) T is the adjoint of x + - 

The condition (2) is also necessary for the existence of spectral factor- 
ization and with a suitable constraint on x + (0) we can ensure that the 
spectral factor x + is unique. 

It is widely known that the solution of several applied problems in 
Control Engineering and Communications require spectral factorization. 
Due to its importance, it is not surprising that, starting from Wiener's 
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efforts [8], a variety of methods has been developed for approximate 
numerical calculation of coefficients 7^ in (4) for a given matrix-function 
(1) (see the survey paper [6]). Most of these algorithms are applicable 
only under the additional restriction on S, e.g., to be scalar or rational. 
There exists an explicit formula for the spectral factor in the scalar case, 

(5) X + {z) = exp jT|^ log S(e*)di) ' 

which is a core of Kolmogorov's exp-log method (see [4], p. 211). There is 
no analog of formula (5) in the matrix case, since in general e A+B 7^ e A e B 
for non- commutative matrices A and B. This is the main reason that 
the algorithms in the matrix case are significantly more difficult. 

An absolutely new approach to the matrix spectral factorization prob- 
lem was proposed by G. Janashia and E. Lagvilava in [3] without im- 
posing on matrix-function (1) any additional restriction apart from the 
necessary and sufficient condition (2). This is the first time that the the- 
ory of Hardy spaces is used for solution of the problem which naturally 
turned out to be very effective since the problem itself is posed in this 
branch of mathematics. 

A sequence of matrix-functions S n is constructed in Janashia-Lagvilava 
algorithm which approximates S in L\ norm, 

||^)-5(z)|| Li -0. 

Then the explicit spectral factorization of S n is performed 

(6) S n (z)=x+(z)(xt(z)r, 
and it is proved that 

(7) \\xt(z)-X + (z)\\ H2 ^0. 

In the present paper, we intend to obtain some qualitative estimation 
of the closeness of Xn to X + which might be useful for practical compu- 
tations of an approximate spectral factor within a given accuracy. We 
achieve this goal under the certain restriction on (1), namely, 

(8) S~\z) e Li(T), 
and estimate 

(9) WxM-x + (z)\\ Hl = \\xM-x + (z)\\ Ll 

from above (see Theorem 1 below). For example, a trigonometric poly- 
nomial matrix-function S(z) = X^feL-w a k zk , ®k are rxr matrix coeffi- 
cients, without zeroes of the determinant on T, det S(z) > 0, \z\ = 1, 
satisfies this condition. It is well-known that the spectral factor x + ( z ) — 
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J2k=o lkZ k * s analytic polynomial of the same order N in this case, and 
the estimation of (9) can be used to obtain the accuracy of approximately 
computed matrix coefficients 7^, k — 0, 1, . . . , N. 

In the present paper we consider only two dimensional matrices, al- 
though, as Prof. Lagvilava informed us, the method we propose can be 
extended to higher dimensional matrices as well. 

It should be mentioned that we do not meet any type of estimation 
in the rate of convergence in any other above mentioned matrix spectral 
factorization algorithms. 

2. Notation. Let D := {z G C : \z\ < 1} and T := dD = {z G C : \z\ = 
1}. L p (T), p > 1, denotes the class of pth integrable complex functions 
with usual norm || ■ \\i p . For / G Li(T), the nth Fourier coefficient of / 
is denoted by c n (f). 

A(D) denotes the class of analytic functions in D. The Hardy space 

H p = j/ G A(D) : ||/|| Hp = (sup |/(re*)|*dt < 00) ' < 00 j 
and the space 

L+(T) := {/ G L p (T) : c n (f) = forn < 0} 

of the boundary values of functions from H p are naturally identified. The 
"+" superscript of a function / + emphasizes that the function belongs 
to L+. f E Hp is called outer, we denote / G H°, if 

f(z) = c-expQ-j^ f±^\ g\f(e u )\dty \c\ = 1, z G D. 
A matrix function U(z), z G T, is called unitary if U(z)U*(z) = I, 

T 

where / is the identity matrix of suitable dimension and U*(z) — U(z) . 

Among the several equivalent norms of a matrix-function we select the 
maximum norm of its entries 

||M(z)|| =max||M^)||, 
and M G L p means that G L p . 

3. The Janashia-Lagvilava algorithm. In order to formulate exactly 
our result, we need to describe the Janashia-Lagvilava spectral factoriza- 
tion method [3] in more details. 



Given a positive definite matrix-function 

'a(z) b(zY 

(10) S(z) = 1 



b(z) c(z) 
having the properties 

(11) a(z) , b(z) , c(z) G Li(T), 

(12) a(z) > 0, a(z)c(z) - \b{z)\ 2 > for a.a. zeT, 

(13) log A(z) := log (a(z)c(z) - |6(^)| 2 ) G L^T), 
the lower-upper triangular factorization is performed at first 



(14) S(z) 



g + {z) \ g+(z) <p(z) 
<p(z) f + (z)J { T r iz) / 



with the scalar spectral factors of a(z) and A(z)/a(z), respectively, on 
the diagonal, i.e. 



(15) g + (z)g+(z) = \g + (z)\ 2 = a(z), 

A(z) 



(16) f + (z)f+(z) = \r(z)f 



a(z) 



and (p(z) = b(z)/g + (z) in (14). The relations (11)-(13) imply that 
loga(z) G Li(T) as well, so that the scalar spectral factorizations in 
(15) and (16) exist. 

The equation |y?(-2)| 2 + \ f + (z)\ 2 = c(z) in (14) and the relations (15) 
and (11) imply that g + , f + and ip are square integrable functions, namely, 

(17) max{||(? + |U 2 , ||/ + || i2 , \\<p\\ L2 } < y^SfZ. 
Thus 

(18) g + ,f + eH°, and^GL 2 (T). 
For certainty, it is assumed that 

(19) g + (0) > and / + (0) > 0. 

Next the function f is approximated in L2 norm by its Fourier series 

00 

<Pn(z) = ^ C k{v)z n , 
k=—n 

so that 

(20) \\ip n {z)-ip{z)\\ L2 ^0, 



and the matrix function S(z) (see (10), (14)) is approximated in L\ norm 
by 



SJz) 



9 + (z) \ /«/+(*) 

<p n {z) f + (z)J { o 



(21) }&(z) 



As it was mentioned in Introduction, an explicit factorization (6) is per- 
formed in [3] by proving that 

o \ / at{z) ti( z y 

^ n (z) f+(z)J \-fflz) 

where the second multiplier in (21) is a unitary matrix-function with 
determinant 1, 

(22) l«M| 2 + |^ + W| 2 = l, |s| = l, 
and are polynomials of order n, 

n n 

a n( Z ) = Yl ^ alld ^( Z ) = Yl ^ ' 
k=0 k=0 

whose coefficients can be found by solving a certain system of linear 
algebraic equations of order n (see [3], the system (17) therein). 
The equations (21) and (22) provide that (see [3], (21) therein) 

(23) det X +(z)=g+(z)f + (z), z e D, n = l,2,... 

Imposing all the same uniqueness conditions on matrices x^(0), n = 
1,2,.. ., say, to be a lower triangular with positive diagonal entries, 

(24) [ x +( )] 12 = 0, k+(0)] u >0, [x+(0)b>0, 
the convergence (7) is proved (see [3], Sect. 3). 

4. The convergence rate estimation. The closeness of the approxi- 
mate spectral factor x+ to x + is estimated in the following 

Theorem I. Let S,S n ,x + ,xt be as in Section 3. Then 

(25) Wxt ~ X + lk < v^ll^ll^ll^-l^ll^ - V \\l 2 + \\<p n - V \\ Ll . 

Proof. We suppress the argument z in what follows. The equation (22) 
implies 

(26) \\a+\\ Loo < 1 and ||/3+|| Loo < 1 n = l,2,.... 
Taking the product in (21), we get 
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The claim of the algorithm that Xn * s a spectral factor implies that 
(28) tt+ and tt+ G L+. 

The fact that and ^ n are square integrable follows also from (18) 
and (26). 

Writing x+ in the similar form as xti m (27) and taking the difference, 
we get 

\\Xn -Xmhi =max{||^ + (a+ -a+)|| Ll , \\g + (P+ - P£)\\ Ll , 
\\(p n (a£ - a+) + (cp„ - (p m )a+ - f + {(3+ - P+)\\ Ll , 

\\<Pn(P: ~ fa) + (<Pn ~ <Pm)P: + f + (< ~ <)|k} 



(30) 



Thus, by virtue of the Cauchy-Schwartz inequality, H/pHli < ||/||l 2 x 
||<?||l 2 , (17), and (26), we have 
(29) 

llXn-XmlUi < VPIk m ax{||«+-tt+|| L2 , ||/?n -/3^||L 2 } + |kn-^mlUi- 

Consider now the second rows of xti an d Xm- 

' <Pn(*+ - f + P+ = , f " /+£+ = *L , 

. + /+<=*+, 1 y? m/ 3+ + /+< = * 2 + m . 

It follows from (30) that 

>„(<*+ - a+) - /+(£+ - £+) = *+ - *+ m - fo> B - y, m )a+ , 

.MPn ~ fa) +f + {<~<) = ^tn ~ ~ (<Pn ~ Vm)fa ■ 

Subtracting the first equation times (/3+ — /3+) from the second equation 
times (a+ — a+), we get 

/ + (|a + - a + m | 2 + \fa - /?+| 2 ) = (*+ - *+J(a+ - «+)- 

(*+ - - /£) - fan - ¥> m )(/£a+ - faO- 

Hence 

(31) |a+ - c*+| 2 + |0+ - 0+| 2 = -!(*+ - *JJ(o+ - «+)- 

- " 0+) " - ^)(/9 + a+ - 0+C+). 

Since 

1 / c — &\ 




(see (10), (8)) and = f (see (16)), we have 

( 32 ) ll^llL<lls _1 IUi<°°- 



Thus, by virtue of (32), (28), (18), (26), and the Cauchy-Schwartz in- 
equality, the summands of the right-hand side expression in the equation 
(31) are integrable. The generalization of Smirnov's theorem (see [5], 
p. 109) claims that if $ = $i/<E>2, where $1 G H 2 and $ 2 G H$ , and 
the boundary values of <3> belongs to Li(T), then $ G ifi. Thus the 
first two functions on the right-hand side of the equation (31) are from 



Hi = Lf(T). Note also that the left-hand side function is positive. We 
can write 

(33) c {K - atf + \fa - faf} = c {^tn ~ n m )K ~ <)}- 

c {^(K t ~ - /£)} - coij+bn - <Pm)(fto£ - /£«+)}. 

Now, it follows from (27), (24) and (19) that 

(34) a+(0) > and /3+(0) = 0, n = l,2,..., 
so that 

(35) c {-^ ((*+ - - fa)} = 

J^((*tn(0) - *£n(0))(#(0) - /3+(0)) = 0. 

Since det(x+(0)) = <? + (0)/+(0) = g+(0)a+(0)V+ n (0) (see, respectively, 
(23) and (27), (34)), we have ^(0) = /+(0)/a+(0), and analogously 
#+ m (0) = /+(0)/a+(0). Thus, taking into account (34), 

(36) c {-^(*t n -*JJ« -«+)} = 

' (^n(0)-^ m (0))(«+(0)-«+(0)) = 



/ + (0) 



1 ^^)-<(o^- ^?:^ )|2 <o. 



«+(0) «+(0)7 v " w mV " «+(0)a+(0) 
It follows from (33), (35), (36) and (26) that 

CoiK - a+f + \fa - fa]*} < 

1 || 1 



| C 0{ — (Vn-<Pm)(Pn<Xn-Pm<Xm)}\ < ~ \ \ ^ ~ 1 1 L , / 



Hence 



l^-«+llL + ll/^-/UHi 2 

'7 1 



27rc {|a+ - a+| 2 + - /3+| 2 } < 2||^|| i2 ||(^ n - ^ 



l 2 



and, taking into account (32), we get 
(37) 



max {||a+ - a+|| i2 , - /3+|| L2 } < (2vfF%J(^ n - y? r , 



2 

l 2 ; 
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It follows now from (29) and (37) that 

\\Xn - XmUi < V^||5 , ||i / 1 2 ||S ,-1 ||J/ 4 |^ n - <Pm\\L? + \\<Pn ~ fmW^ 

and by virtue of (7) and (20), we get (25). □ 
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